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Abstract 

We show that, in the framework of Caratheodory's approach to thermodynamics, 
one can implement black hole thermodynamics by realizing that there exixts a quasi- 
homogeneity symmetry of the Pfafnan form 5Q rev representing the infinitesimal heat 
exchanged reversibly by a Kerr-Newman black hole; this allow us to calculate read- 
ily an integrating factor, and, as a consequence, a foliation of the thermodynamic 
manifold can be recovered. 
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1 Introduction 



We consider black hole thermodynamics in the framework of Caratheodory's 
approach to thermodynamics, which postulates the integrability of the Pfaf- 
fian form d~Q rev representing the infinitesimal heat exchanged reversibly [1,2]. 
The integrability of SQ rev means that there exist an integrating factor /i and 
a function a (called "empirical entropy" [1]) such that 5Q rev = fida. In stan- 
dard thermodynamics, the integrability of 5Q rev is a consequence of Clausius 
inequality, which ensures that the absolute temperature T is an integrating 
factor for 5Q rev , in such a way that 5Q rev = TdS, where 5* is the entropy 
("metrical entropy" [1]). The existence of the entropy function is part of the 
second law (the principle of entropy increase for thermodynamic processes of 
closed and thermally insulated systems is the other part [2]; it is related to 
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the strict inequality, which holds in irreversible processes, in Clausius inequal- 
ity). If the integrable Pfaffian form 5Q rev displays a symmetry, in a sense to 
be described in the following, it is also possible to calculate explicitly and 
readily an integrating factor by means of elementary tools of differential ge- 
ometry [3,4,5]. As a consequence, an explicit construction of the foliation of 
the thermodynamic manifold into disconnected adiabatic hypersurfaces satis- 
fying 5Q rev = is readily calculated. 

By postulating a natural form for 5Q rev in the case of black holes of the Kerr- 
Newman family, we can introduce a notion of temperature and of entropy 
for black holes without referring a priori to the laws of black hole mechanics, 
thanks to the integrability of 5Q rev and to the presence of a quasi-homogeneity 
symmetry of the Pfaffian form SQ rev . Particularly, we can generate a poten- 
tial which is then related to the entropy of the black hole. Both the entropy 
and the temperature appear as derived quantities. We point out also that the 
following construction for black holes, being made an explicit use both of the 
integrability condition and of the symmetry of SQ rev , is of interest also for 
thermodynamicists. 



2 Pfaffian form and symmetry 

We recall that, for a Kerr-Newman black hole, due to the no-hair theorem, 
the only parameters available are M, Q, J (mass, charge and angular momen- 
tum); for an insulated black hole they play the role of conserved charges, and 
are chosen as independent variables in the thermodynamic domain, which is 
assumed to be the non-extremal manifold M 4 — M 2 Q 2 — J 2 > (the extremal 
sub-manifold M 4 — M 2 Q 2 — J 2 = is a boundary of the former, and is tem- 
porarily not taken into account. Some more discussion on this topic is found 
in sect. 4). We look for a natural infinitesimal reversible form for the first 
law. Then, we have to find out a Pfaffian form such that the first law in its in- 
finitesimal form is implemented, and moreover, such that it is integrable (thus, 
the first part of the second law is ensured). In the case of a charged rotating 
system, it would be natural to define SQ rev = dU + pdV — (fidq — ujdl — /idN, 
where q is the electric charge, / is the angular momentum, <fi is the electrostatic 
potential and ui is the angular velocity. We introduce then the angular velocity 



n = 



j i 



(i) 



M 2M 2 -Q 2 + 2M^M 2 - Q 2 - J 2 /M 2 



and the electric potential 




(2) 
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of the black hole. Both f2 and $ can be assigned on a purely geometrical foot- 
ing, without any a priori knowledge of black hole thermodynamics. Analogy 
with standard thermodynamics leads us to associate infinitesimal variations 
of Q, J with work terms — $g?Q and —VtdJ; infinitesimal variations of M can 
naturally play the role of dU in standard thermodynamics, where U is the 
internal energy. We then define 

5Q rev = dM- <5>dQ - ttdJ. (3) 

Definition (3) is natural, in fact the (rest) mass can be identified with (a term 
of) the internal energy (the rest mass of a fluid can be considered as a term 
of the internal energy in standard thermodynamics; see e.g. [6]); moreover, 
as seen, the work terms appear as standard work terms. The most evident 
difference with the case of a standard system consists in the absence, for the 
black hole case, of the pdV term. This lack is associated with the lack of a 
notion of volume in the black hole well known. Even a pressure cannot 

be defined, and the same is true for the particle number N. 
It is stressed that the infinitesimal variation dM — QdQ — VtdJ is taken along 
stationary black hole solutions of the Kerr- Newman family, because of (1), (2); 
this means that the Einstein equations are satisfied for each state involved 
in the aforementioned variation. Moreover, these solutions of the Einstein 
equations are considered as black hole equilibrium states, to be compared 
with equilibrium states of standard thermodynamics. 

The Pfaffian form SQ rev is everywhere non-singular, i.e., there is no point of 
the thermodynamic domain where all the coefficients of the differential form 
vanish (this property holds also for points of the extremal boundary). It is easy 
to show that 5Q rev is smooth on the non-extremal manifold and is completely 
integrable, that is, it satisfies the condition 5Q rev A d(5Q rev ) = 0, i.e. 

- dj$ + d Q n + <s>d M n - nd M ® = o. (4) 

Being 5Q rev a one-form in three variables, this integrability condition is surely 
non-trivial (it would be trivial in the case of two variables). Notice that a 
different choice for the sign of the work terms in (3) would lead to a non- 
integrable Pfafnan form, as it is easy to verify. In the appendix, we extend 
the present considerations to the case where a magnetic monopole charge P is 
allowed. It is remarkable that integrability is not postulated, but it simply fol- 
lows from considering infinitesimal "on shell" variations, i.e. variations along 
the aforementioned solutions. 

We can also find an integrating factor by using the quasi-homogeneity symme- 
try of the Pfaffian form (3). 1 In fact, under the quasi-homogeneous transfor- 

1 A Pfaffian form w = Ya=i u i( xl -> • • • > x n )dx % is quasi-homogeneous of degree r G IR 
and weights (a±, ... ,a n ) € IR™ if, under the scaling x 1 , . . . , x n i-> A^x 1 , . . . , \ an x n 
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mation [7] (also called "similarity transformation" and "stretching transfor- 
mation" ) M i-> \ a M; Q ^ \ a Q; J ^ \ 2a J, one obtains 5Q rev i-> \ a 5Q rev , 
i.e., 5Q rev is quasi-homogeneous of degree a. (a, a, 2a) are defined to be the 
weights of M,Q, J respectively and they have to be determined. Let us de- 
fine the so-called Euler vector field [7], which is infinitesimal generator of the 
transformation 

D ^ aM JL + aQ JL + 2aJ jL, (5) 



we have introduced above a label a which underlines that a is not yet fixed 
{{Da} a is a one-parameter family of Euler vector fields); let the correspond- 
ing Lie derivative be £_d q ; then, the quasi- homogeneous transformation is a 
symmetry for 5Q rev (see e.g. [8,9]), in the sense that 

(L Da 5Q rev ) A 5Q rev = 0. (6) 



In fact, L Da 5Q rev = a8Q rev . An integrating factor f a such that the form 
8Qrev/fa is exact is f a = iD^Qrev — 8Qrev(D a ). For a proof that f a is 
an integrating factor see Ref. [3] and, for the homogeneous case, see e.g. 
Ref. [9] and also Ref. [4], where an application to ordinary thermodynam- 
ics can be found. In our case, one obtains f a = a(M — $Q — 2QJ), explicitly 
fa = ct\J M 2 — Q 2 — J 2 /M 2 , which is not identically vanishing, thus D a is 
associated with a transversal (or non-trivial) symmetry [3] (i.e., D a does not 
belong to the distribution of codimension one associated with the kernel of 
SQrev)- We remark that the integrating factor f a is proportional to the hori- 
zon coordinate c introduced by B.Carter in [10]. Then, f a is constant on the 
horizon. 



3 Foliation of the thermodynamic manifold 



Frobenius theorem for the Pfafnan form SQ rev on the non-extremal manifold 
can be invoked and a foliation of the non-extremal manifold can be generated 
thanks to the integrability property (4). The non-extremal thermodynamic 
space is foliated by the submanifolds (of codimension one) which are solutions 
of the Pfafnan equation 5Q rev = 0. The leaves of the foliation of codimension 
one are surfaces where the potential associated with SQ rev /f a is constant. Let 
(Mo, Qo, Jq) be a reference state and T be any path connecting the reference 

one finds u> i— > \ t oj. This happens if and only if to^x 1 , . . . , x n ) is quasi-homogeneous 
of degree = r — ati for alH = 1, . . . , n, which means that Ui(\ ai x 1 , . . . , X an x n ) = 
\ r - a *Wi(x 1 ,...,x n ). 
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state to the state (M, Q, J) of interest. By choosing e.g. a rectangular path 
(M , Qq, J ) -> (M, Q , J ) -> (M, Q, J ) -> (M, Q, J) contained in the non- 
extremal manifold, one finds 



S a (M,Q,J)-S a (M ,Q ,J ) 




M 2 b 2 (M,Q,J) + J 2 /M 2 \ 
M 2 b 2 (M ,Q ,J ) + J 2 /M 2 ) 1 } 



where b(M, Q, J) = (1 + yjl - Q 2 /M 2 - J 2 /M 4 ). The argument of the loga- 
rithm is proportional to the black hole area A = 4ir(M 2 b 2 (M, Q, J) + J 2 /M 2 ). 
We have generated a foliation of the parameter space of Kerr-Newman black 
holes. The leaves are the surfaces A = const., as expected, but we cannot yet 
determine the so-called metrical entropy [1] in the case of black holes. 
We now introduce an assumption which requires some discussion. The above 
procedure is a generalization, discussed in Ref. [3], of the procedure one can 
develop for standard thermodynamics [4]. In the case of standard thermody- 
namics of homogeneous systems, the Pfaffian form 5Q rev = dU + pdV — fidN 
in Gibbsian variables (U, V, N) is homogeneous (for the definition of homo- 
geneous differential form see e.g. Ref. [9]). The generator of the symmetry is 
the "Liouville" operator Y = Ud v + Vd v + Nd N and the integrating factor is 
SQ rev (Y) = U + pV — jiN . For a standard thermodynamic system one finds 
that dS = SQ rev /f = dS/S, where S is an extensive function which coincides 
with the metrical entropy of the system and corresponds to the fundamental 
relation in the entropy representation [4]. This deduction is corroborated by 
appealing to the homogeneity of S in Gibbs' approach, which allows to find 
TS = U + pV — fiN = 5Q rev (Y), i.e., the integrating factor coincides with TS. 
We proceed by analogy with the formalism of thermodynamics just sketched, 
which means that we assume that the metrical entropy is the unique 2 quasi- 
homogeneous function S of degree one which satisfies dS = SQ rev / f = dS/S. 
We refer to [3] for a proof that such an S exists and is unique. This assumption 
about the role of S such that 5Q rev / f = dS/S is correct both in the case of 
standard thermodynamics in the Gibbs space [4] and in the case of black hole 
thermodynamics, and the latter case is analyzed in the following. 
The potential S a such that dS a = dS a / S a for the black hole case is 

S a = c a A^ 2a (8) 



where c a is an undetermined constant. We have a one-parameter family of 
possible metrical entropies (and also fundamental relations in the entropy 
representation) which satisfy D a S a = S a , in analogy with YS — S of standard 

2 Uniqueness holds within a multiplicative constant. See [3]. 



5 



thermodynamics. Our result (8) agrees with the result contained in Ref. [11] 
but we work in a more general framework where no reference to the laws 
of black hole mechanics is made [notice also that in our expression for S 
no additive constant appears, due to quasi-homogeneity symmetry]. There 
is still an ambiguity due to the undetermined value of a, which means that 
we know the ratio between the weights of M, Q, J but not yet the weights 
themselves. Notice that this ambiguity does not occur in the case of standard 
thermodynamics, where the weights of (U, V, N) are known and they are all 
equal to one. We recall that the metrical entropy is assumed to belong to the 
one-parameter family {S a } a . For each a the temperature is T a = (dS a / dM)^ 1 
and it is a quasi- homogeneous function of degree a — 1 and weights (a, a, 2a). 
It is useful to realize that 

s ^(^ (c '^ )1/2 ^(^ (Sl ' 2)1/2 " (9) 

r a . = 2 a (ci/2> ' /2 ° (s 1/2 ) 1 - 1 / 2 T 1/2 . (10) 

For any a one gets T a dS a = dM — §dQ — QdJ. The black hole area is known 
to be a superadditive function of M, Q, J. Superadditivity of the entropy, 
which plays a fundamental role when one considers the merging of two black 
holes, does not fix a. 3 Cf. also [11] for a discussion concerning an analogous 
ambiguity (our I /2a is 7 therein). The Hawking effect is necessary in order 
to give us an actual thermodynamic meaning to our calculation; it also fixes 
a, in fact, in order to identify the temperature of the black hole with the 
Hawking one it is mandatory to choose a — 1/2. There is also a multiplicative 
constant (namely, Ci/ 2 ) which has to be determined. By comparison with the 
Hawking effect, one finds that C1/2 = 1/4. The above ambiguity can also be 
resolved phenomenologically, in perfect agreement with the phenomenological 
nature of a thermodynamic approach [3]; one should determine M,Q, J and 
then plot T(M,Q, J) from measurements of the temperature, a = 1/2 (and 
also C1/2 = 1/4) should come out again. 

It is remarkable that, as a consequence of the quasi-homogeneity of black hole 
entropy, one gets the following generalized Gibbs-Duhem equation, which is 
analogous to the Gibbs-Duhem equation of standard thermodynamics: 

M 2 d ( -M - Q 2 d ( -%=) - Jd (%) = 0. (11) 



3 0<a<l/2isa sufficient condition for preserving superadditivity. 
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This follows from S = ioi^Qrev/T) and from dS = 5Q rev /T, where D = Di/ 2 - 
It is easy to show that this generalized Gibbs-Duhem equation corresponds to 

- iD d(*^)=0. (12) 



In fact, dS = -i D d(5Q rev /T) + L D (5Q rev /T) = -i D d(5Q rev /T) + 5Q rev /T; 
from the latter equality and from dS = 5Q rev /T equation (12) follows. A sim- 
ple rearrangement of the terms one obtains by making explicit (12) gives then 
(11). See [3] for a general setting and [12] for the case of standard thermody- 
namics. 

Contrarily to the naive expectation, the laws of black hole mechanics give 
no unique hints about the value of a, they don't fix uniquely the metrical 
entropy and the absolute temperature of the black hole. For any a one gets 
T a dS a = dM — QdQ — QdJ, to be compared with the differential form of 
the first law. Moreover, one finds that f a = T a S a , which implies a(M — 
$Q — 20, J) = T a S a = 2aTi/ 2 Si/2- By comparison with the first law in the 
finite form one realizes that Ty 2 = fc/(87rci/ 2 ). The choice of a generic a is 
equivalent to the the substitutions A i— > A a and k i— > k a , where A a = A x l 2a 
and k a = 2ak/A 1 / 2a ~ 1 , which implement both the differential form and the 
finite form of the first law (the latter appears as k a A a = 8ira(M — $Q — 2QJ) 
which is equivalent to the well-known one). Notice that k a is constant on the 
horizon, thus the zeroth law of black hole mechanics is not sufficient in order 
to select a — 1/2. 



4 The extremal boundary 



The extremal submanifold is very problematic. It is easy to show that 5Q rev = 
on the extremal submanifold, i.e. the extremal submanifold is still an inte- 
gral submanifold of the Pfaffian form [13]. Nevertheless, there is an important 
property which fails in the case of states belonging to the extremal submani- 
fold. In fact, given a point of the extremal submanifold, there exist two kinds 
of adiabatic paths having the given state as initial point. One is a path lying 
on the extremal submanifold, the other is an "isoareal" path, i.e. a path start- 
ing from the extremal submanifold and reaching non-extremal states each of 
which has the same area as the initial extremal state [13]. In absence of the lat- 
ter class of solutions, the extremal states would represent a leaf of a foliation, 
thus they would be adiabatically disconnected from the non-extremal states. 
Instead, their presence can jeopardize the second law of thermodynamics. A 
detailed discussion of this topic and of the third law in black hole thermo- 
dynamics is the subject of Ref. [13]. See also [14] for the case of standard 
thermodynamics. 
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5 Conclusions 



The approach to black hole thermodynamics by means of Pfaffian forms we 
have discussed (Caratheodory's formalism) represents a further corroboration 
of the fact that black hole thermodynamics is a form of thermodynamics, even 
if to large extent exceptional. Quasi-homogeneity symmetry of SQ rev plays an 
important role in allowing to calculate an integrating factor and, then, to gen- 
erate a thermodynamic potential depending on a undetermined parameter a 
which can nevertheless be fixed phenomenologically, as seen. Notice that our 
approach can be extended in a straightforward way to KN-AdS black holes 
[15]. Also in this case, there is a quasi-homogeneity structure in the Pfaffian 
form, as it can be easily realized. It is remarkable that quasi-homogeneity is a 
symmetry occurring also in the thermodynamics of other self-gravitating sys- 
tems, like non-relativistic fermionic matter and self-gravitating radiation. In- 
dependent thermodynamic variables and their weights change 4 , but a quasi- 
homogeneous symmetry appears again. See [3] on the latter topic. The corre- 
sponding lack of homogeneity can be related to the purely attractive nature 
of gravity. 

This kind of thermodynamic approach can be insightful also from the point 
of view of a more general discussion concerning horizon thermodynamics. The 
availability of a meaningful 5Q rev allows to discriminate between cases where 
there is a complete thermodynamic structure at hand, which can be associ- 
ated with the Einstein equations (e.g. because 5Q rev is integrable "on shell") 
and cases where Hawking temperature seems to be related simply to kine- 
matics. Cf. [16]. The case of black holes belongs to the former class, a full 
thermodynamic structure exists and the role of the equations of General Rel- 
ativity in ensuring the laws of thermodynamics, enhanced e.g. in Ref. [16], is 
corroborated in this framework. Notice that, from this point of view, the be- 
havior of General Relativity is, to some extent, intermediate with respect to a 
macroscopic phenomenological approach, like classical thermodynamics, and 
a microscopic approach, like statistical mechanics. A macroscopic ("thermo- 
dynamic") point of view is adopted in treating variables like M,Q, J; on the 
other hand, field equations furnish $, Q and ensure an integrability condition 
which, for standard systems, should be an outcome of statistical mechanics 
[statistical mechanics should calculate the analytic form of the functions Q 
(they are phenomenological interpolations for thermodynamics); moreover, it 
should justify an integrability condition which is only a postulate in stan- 
dard thermodynamics; statistical mechanics should allow to determine both 
the metrical entropy and the weights of the variables M, Q, J [3]]. As far as a 
"cosmological horizon" like De Sitter (one parameter) is concerned, it corre- 



4 One has (U, V, N) as independent variables in the fermionic matter case, and 
(U, V) in the self-gravitating radiation case. 
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sponds to a solution of General Relativity and 5Q rev can be given (but a too 
tight thermodynamic space does not allow a non-trivial integrability condi- 
tion). On the other hand, in the "acceleration horizon" case (like e.g. Rindler 
case) there is a too poor thermodynamic structure, in the sense that there is 
no first law. The same comment holds true in the case of "acoustic horizons" 
(no first law is known). 



A Magnetically charged black holes 

An extension is represented by a rotating charged black hole: 

SQ rev = dM -QdQ-VdP-ndJ, (A.l) 



where P is the magnetic monopole charge and 



^=T7 r ( A - 2 ) 

M 2M 2 -Q 2 -P 2 + 2M^/M 2 -Q 2 -P 2 - J 2 /M 2 

Q(M + JM 2 - Q 2 - P 2 - J 2 /M 2 ) 

$ = 1 , (A.3) 

2M 2 - Q 2 - P 2 + 2M^M 2 -Q 2 -P 2 - J 2 /M 2 



P{M + JM 2 -Q 2 -P 2 - J 2 IM 2 ) 

* = . =. (A.4) 

2M 2 - Q 2 - P 2 + 2MyJM 2 - Q 2 - P 2 - J 2 /M 2 

Notice that 5Q rev is "on shell", i.e., it is taken along (stationary) solutions 
of General Relativity. 5Q rev A d(5Q rev ) = corresponds to four integrability 
conditions: 



l MQP = -dp® + d Q ^ + <$>d M ^ - ^d M <$> = 0, (A.5) 

l MQJ = -dj$ + d Q Q + <5>d M Q - fid M $ = 0, (A.6) 

i MPJ = -dj^ + d P n + vo M n - nd M v = o, (A.7) 

l QPJ = $(<9j^ - d P Q) - ^(,9/$ - d Q Q) + n(d P $ - d Q V) = 0. (A.8) 

Even in this case, the integrability is verified and SQ rev is a quasi-homogeneous 
Pfafnan form. The Euler vector field is given by 

1 d 1 d 1 d d 

2 dM 2^dQ 2 dP dJ K J 



The integrating factor is easily found, and, again the area law can be obtained. 
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